Generation of spin currents and spin densities in systems with reduced symmetry 
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We show that the spin-current response of a semiconductor crystal to an external electric field 
is considerably more complex than previously assumed. While in systems of high symmetry only 
the spin-Hall components are allowed, in systems of lower symmetry other non-spin-Hall compo- 
nents may be present. We argue that, when spin-orbit interactions are present only in the band 
structure, the distinction between intrinsic and extrinsic contributions to the spin current is not 
useful. We show that the generation of spin currents and that of spin densities in an electric field 
are closely related, and that our general theory provides a systematic way to distinguish between 
them in experiment. We discuss also the meaning of vertex corrections in systems with spin-orbit 
interactions. 



Ground-breaking work in past years has turned semi- 
conductor spin electronics into a richly rewarding field, 
both theoretically and experimentally. The application 
of an electric field to a semiconductor sample gives rise 
to a nonequilibrium spin current P, 0, |j| as well as 
a steady-state spin density in the bulk of the sample 

[fi 2_ 3 • The latter was first observed several decades 
ago 5|, while the recent imaging 0] and direct measure- 
ments of spin currents 0, [13] , together with the achieve- 
ment of the room-temperature spin-Hall effect ll| have 
brought the promise of spin-based electrical devices closer 
to fruition and stimulated an enormous amount of re- 
search Other schemes for generating and detecting 
spin currents have been implemented or proposed (l3| . 

Aside from the specter of technological advances, the 
development of electrical spin manipulation techniques 
has brought the fundamental ph ysic s of spin transport 
under intense scrutiny [li, [11, [11 [3, [l^, [211 De- 

bate has focused on definitions of spin currents [14| , on 
whether spin currents are transport or background cur- 
rents , whether scattering-independent or scattering- 
dependent contributions are dominant [l3 |. on the role 
played by spin Coulomb drag [l9| . and on spin accumu- 
lation at the boundary [sO]. Whereas most studies to 
date have paid significant attention to common semicon- 
ductors and asymmetric quantum wells 22, [23|, [2J, |25|, 
H, [H, [H, [H [aa, [m, [H [si, [H [ail , recent developments 
call for an in-depth investigation of neglected aspects. 

In this article, we first discuss the relationship between 
spin currents in a crystal and the symmetry of the under- 
lying lattice. Such an analysis has been enlightening in 
the context of nonequilibrium spin densities excited by an 
electric field. If the response of the spin density s to an 
electric field E is given by s" = Q'jEj, nonzero compo- 
nents for the material-specific spin density response ten- 
sor QJ are permitted only in gyrotropic crystals ^ . Yet 
such an analysis has not been done for spin transport. We 
therefore determine the components of the spin-current 
response tensor allowed by symmetry in an electric field, 
providing systematic proof that spin currents in response 



to an electric field can be much more complex than the 
spin-Hall effect [s^ . This result is completely general and 
is not sensitive to the definition of the spin current or to 
whether the electric field is constant or time-dependent. 
The subsequent calculations confirm these predictions by 
determining the corrections to the density matrix present 
in an electric field and demonstrating their intimate re- 
lationship with spin precession. We argue that, if spin- 
orbit interactions are present only in the band structure, 
only one contribution to the spin current exists, which 
appears intrinsic in the weak momentum-scattering limit 
and extrinsic in the strong momentum-scattering limit. 
We show that spin currents and bulk spin densities in 
an electric field arise from linearly independent contribu- 
tions to the density matrix, and that certain setups can 
distinguish between them and measure effects due solely 
to spin currents. Our work considers realistic scattering 
potentials and is very relevant to experiment, where the 
high symmetry often assumed in theoretical approaches 
is usually not present. 

The spin current operator is defined as = ^{s"'vi + 
iiis'^), where s" represents the spin component a, and the 
velocity operator Vi — l/h{dH/dki), with H the Hamil- 
tonian. An alternative, more realistic definition of the 
spin current has been proposed [l^ . according to which 
— d/dt {fis'^). Yet we remark that from a symmetry 
point of view these two definitions are equivalent so that 
the following symmetry analysis applies to both defini- 
tions. The spin current JT" is a second rank tensor that 
can be decomposed into a pseudoscalar part, an anti- 
symmetric (spin-Hall) part, and a symmetric part. The 
pseudoscalar part is tr(j7') = ^ s-v and represents a spin 
flowing in the direction in which it is oriented. The sym- 
metric and antisymmetric parts are given, respectively, 
by i(s'^Wi±'0(jS*). The pseudoscalar and symmetric parts 
will be referred to as non-spin-Hall currents. Under the 
full orthogonal group only the antisymmetric (spin-Hall) 
components of JT" are allowed, indicating that these com- 
ponents are always permitted by symmetry. 

In general, the spin current response of a crystal to an 
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electric field E is characterized by a material tensor T 
defined by J°' = T^^Ej . For the 32 crystallographic point 
groups the symmetry analysis for the tensor T is 
established by means of standard compatibility relations 
(sst - One is particularly interested in those groups in 
which non-spin-Hall components may be present. Our 
calculations show that the pseudoscalar part of the spin 
current is only allowed by 13 point groups, while the 
symmetric part is allowed in all systems except those 
with point groups O, Td (zinc blende), or Oh (diamond). 
Lower symmetries, allowing non-spin-Hall currents, are 
characteristic of systems of reduced dimensionality. 

It is well known that the generation of a spin density 
by an electric field is the inverse of the circular photo- 
galvanic effect [1, [13] , while the spin-Hall effect also has 
an inverse [T^]. Since spin densities induced by electric 
fields are restricted to gyrotropic crystals \^ the same 
restriction applies to the circular photogalvanic effect. 
Similarly, the symmetry analysis developed here is also 
applicable to the inverse spin-Hall effect which needs to 
be complemented by an inverse non-spin-Hall effect in 
systems with reduced symmetry. 

In order to verify the proposition that non-spin-Hall 
currents may exist in a variety of crystals, we discuss spin 
currents induced by an electric field in spin- 1/2 electron 
systems. For electrons the effective Hamiltonian is writ- 
ten as H ^ i?kin + Hso, where i/kin = ^^fc^/2m* and 
^^so = 5 ■ where SI is a momentum-dependent ef- 
fective Zeeman field. In the weak momentum-scattering 
regime we have epTp/h ^ D,Tp/h ^ 1, with Tp 
the momentum relaxation time and ep the Fermi en- 
ergy, whereas in the strong momentum-scattering regime 
SFTp/h:$> 1 > QTp/h. 

The system is described by a density operator p, which 
is expanded in a basis of definite wave vector as 



P = 
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A constant uniform electric field E is included in the 
crystal momentum through the vector potential A such 
that k — q + eA/h, and the wave functions are chosen 
to have the form \tpnk{t)) = e'l"^ \unk{t)) , where \unk{t)) 
are lattice-periodic functions that are not assumed to be 
eigenfunctions of the crystal Hamiltonian. In this basis 
the matrix elements Pnn'kk' (t) form the density matrix 
and the impurity potential has matrix elements Ukk' 1 + 
Vfcfc', where Vkk' is the spin-dependent part. 

The time evolution of the density operator is given 
by the quantum Liouville equation, which allows us to 
derive rigorously the time evolution of the part of the 
density matrix p = Pnn'kk diagonal in wave vector. We 
subdivide p = po + PEi where po is given by the Fermi- 
Dirac distribution, and the correction pE is due to the 
electric field E. To first order in E, pE satisfies 

dpE i ,rj ] , 9r \ eE dpo , . 

-g^ + j:[H.pE] + J{pE)^j--^, (2) 



where J (pe) is the collision integral to be given below. 
The eigenvalues of H are e± = Eq ± fl/2, with eo = 
h?k'^/2'm*. The spin current operator is simply — 
tikis'^ /m* + jhOff^ /dki. In electron systems we usually 
have Hso < ^^kin, thus po ~ /o 1 + -^so dfo/deo, where 
/o(eo) is the Fermi-Dirac distribution. pE is divided into 
a scalar part and a spin-dependent part, pE — /e1 + Se- 
To first order in iJso/^^kin the scattering term can be 
expressed as J (pe) = (Jo + Js + Jv) (/s) + Jo (Se), where 



Jo Ue) 



27rn,- 



dfik' 
d'^k' 



JsUe) = ^<T 



^\Ukk'?{fE- fE)5{eo-e'o)ii^ 
\Kkk'\' ifE - /'e) 
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x{{n + n') [6{e+ - eV) - <5(e_ - e'_)] 

+ (0 - Cl') [5{e+ eL) - 5{e^ e'^)]}, (3b) 



Jv Ue) 



27rn,- 



dfik' 



■Ykk'{fE-fE)6{eo-e'o). (3c) 



In the above rii is the impurity density, d is the dimen- 
sionality of the system, primed quantities denote func- 
tions of fc', $7 is a unit vector along $7, and Yfcfc' = 
/ d'^k" /{2TTY{Ukk"Vk"k- + Vkk-Uk"k')- Equation ^ 



is the usual scalar scattering term, Eq. pb)l is due to 
band structure spin-orbit coupling, and Eq. ((3c|) is due 
to spin-orbit coupling in the impurities. We will assume 
henceforth that band structure spin-orbit interactions are 
much stronger than those due to impurities. 

For Je we obtain the standard correction = 
(e E Tp/h)dfo/dk 1. As a result, the effective source term 
that enters the equation for Se is Ss — Js (/b), where Ss 
is the spin-dependent part of {eE/fijOpo/dk. In analogy 
with the customary Gram-Schmidt orthogonalization of 
vectors, this effective source term can be divided into two 
parts, Eg — Js {/e) = '^\\ + of which E|| commutes 
with the spin-orbit Hamiltonian 

tr{[S],-J,(fe)]g,o} 

(4) 

while Ti±_ is the remainder. In matrix language Ex is 
orthogonal to the Hamiltonian, thus tr(Ex Hso) = 0. To 
find E|| and E_l we define projectors Py and P± onto 
and orthogonal to Hso, respectively. Acting on the basis 
matrices <t, Pycr = 2n Hgo/^"^, while P±(Tx = [{^l + 
il^) ax—^x^y o'y—flx^z Cz]/^fej and the remaining terms 
are obtained by cyclic permutations. 

Se can likewise be divided into two linearly indepen- 
dent contributions, S\\ and S±, the former of which com- 
mutes with Hso while the latter is orthogonal to it. It 
is helpful to think of S\\ as the distribution of conserved 
spins parallel to SI and of S± as the distribution of pre- 
cessing spins perpendicular to it. S\\ satisfies 



dS\\ 
'dt 



P\\[Jo{Se)] 



^11 



(5) 



3 



This equation can be solved iteratively for any scatter- 
ing. To solve the equation for S±, on the other hand, 
one needs to expand S± in the strength of the scattering 
potential \U\'^, as 5*^ = S^o + S^i + 0(1^^1"). (It is eas- 
ily shown that the first term in the expansion must be 
zeroth order in while starts at order —1.) The 
equations for these contributions are 

dS^o ■ 



dt 
dt 



[Hso, S±o] 



— [Hso, S±i] 



P^[Jo(5||)], (6a) 



(6b) 



It follows that if Y,± — Pi [Jo {S\\)] vanishes in Eq. ([6a| . 
as it does for the linear Rashba model, then all the con- 
tributions to S± vanish. In general, the equation for S±o 
is also solved iteratively for any scattering, but a closed- 
form solution for Se is not always possible. An enlighten- 
ing closed-form solution is, however, possible for short- 
range impurities, where Jo{Se) = {Se — Se)/tp, with 
the bar denoting averaging over directions in fc-space and 
Tp = ?i'^/(m* |Z//p). The series of equations of increasing 
order in \U\'^ give two geometric progressions that sum 
to 



'II 



^11 (l-^ll 



Si 



S±rp + Pj_5|| n X {Y,j_Tp + P^S\\) ■ (TTp 



(7a) 



(7b) 



In the general case considered below the complex expres- 
sions for Se, 5] II and S_l will not be given explicitly. 

Our analysis clarifies the relation between steady-state 
spin currents and spin densities that appear when an elec- 
tric field is applied to a semiconductor. Since the spin op- 
erator is even in k and the spin current operator is odd in 
fc, it emerges, after evaluating and ^j,, that the RHS 
of Eq. ((7a|) is responsible for steady-state spin densities 
B 0, [3|, while the second term on the RHS of Eq. ([70) 
gives rise to spin currents P, [3, H, 0] ■ [The first term on 
the RHS of Eq. (j7b[) vanishes in both the weak and the 
strong momentum scattering limits.] We conclude that 
spin densities arise from S\\, the distribution of conserved 
spin while spin currents arise from S±_ , the distribution of 
precessing spin. Thus nonequilibrium spin currents are 
due to spin precession (as first demonstrated in Q) while 
nonequilibrium spin densities 0, |3, are due to the ab- 
sence of spin precession. The physical picture for the 
latter mechanism is as follows. Each spin on the Fermi 
surface precesses about an effective field 17(fc) and the 
spin component parallel to 0,(k) is preserved. In equi- 
librium the average of the conserved components is zero, 
but when an electric field is applied the Fermi surface 
is shifted and the average of the conserved spin compo- 
nents may be nonzero. This intuitive physical argument, 
to our knowledge absent from the literature, explains why 
the nonequilibrium spin density oc and requires seat- 



that, on the other hand, spin currents contain only terms 
ocr^" with n = 0,1,2,.... 

It can be seen from Eq. (j7bp that there is only one 
contribution to the spin current, which in the weak 
momentum-scattering (intrinsic) limit is independent of 
Tp and in the strong momentum-scattering (extrinsic) 
limit is oc Tp. It was also noted that, if the RHS of 
Eq. ([6a| vanishes, then ^j, vanishes to all orders. We 
conclude that, if spin-orbit interactions exist only in the 
band structure, the distinction between intrinsic and ex- 
trinsic spin contributions to the spin current is not useful. 

In calculations of spin currents based on the Green's 
functions formalism vertex corrections play an important 



role 
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28|. The above analysis suggests that con- 
tributions due to two processes are contained in these 
vertex corrections. Scattering renormalizes the driving 
term for 5'^; to — J, (/b), and it mixes the conserved 
and precessing spin distributions, as in Eq. (|Sa)) . This im- 
plies that vertex corrections to S contain the influence 
of a steady-state spin density on J' . As this spin density 
occurs only in gyrotropic materials [H], we expect vertex 
corrections to spin currents to vanish in non-gyrotropic 
materials. 

Engel et al. 31] showed that, when band structure 
spin-orbit interactions are negligible, spin currents arise 
from skew scattering. A comparison of our results with 
those of Ref. [3l[ shows that scattering can give rise to 
spin currents of qualitatively different forms depending 
on whether spin-orbit interactions are present in the band 
structure or not. In the absence of band structure spin- 
orbit interactions, scattering processes involve scattering 
of conserved spins, but in the presence of band structure 
spin-orbit interactions, scattering processes involve the 
scattering of precessing spins. A general analysis of band 
structure spin-orbit interactions and skew scattering on 
the same footing remains to be undertaken. 

For known cases, our theory agrees with previous work. 
2D Hamiltonians with spin-orbit coupling linear in k 
give zer o sp in current for any scalar scatterin g p oten- 



tial [28|, l39|, including short-range [20|, l26|, |27|, l28|, [2 



and small-angle scattering [20|, |2l[. For Hamiltonians 
characterized by one Fourier component N [2lj the spin 
current oc N . In 3D the correction to the spin current 
due to Js (/b) vanishes for the fc'^-Dresselhaus model and 
short-range impurities. Our results also agree with pre- 
vious calculations of spin generation Q. 

As a specific example, we investigate the spin cur- 
rent response tensor T in systems of low symmetry. A 
strong justification for this choice is that experiment of- 
ten studies low-symmetry structures whereas theory is 
often done for high-symmetry models. To reduce the 
symmetry of a system one can either apply strain or 
lower its dimensionality. In strained bulk samples, non- 
spin-Hall currents are expected to be proportional to the 
strain tensor and therefore small. We concentrate thus 



tering to balance the drift of the Fermi surface. We note on quantum wells and, considering for definiteness a sym- 
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metric 150-A wide GaAs well grown along [113]. The 
axes of the coordinate system are x = [332], y = [110] 
and z ~ [113]. Spin-orbit interactions are described 
by A:-linear and /c'^-Dresselhaus terms (40| and impurity 
scattering by a screened Coulomb potential, with the 
Thomas- Fermi wave vector — Akp at carrier density 
n = 1.6 X 10^^ cm~^. We do not consider skew scattering, 
which is important when band structure spin-orbit inter- 
actions are weak 11, 3lj. In units of e/(87r), we obtain 
= 0.398, T^y = 0.12, T^^ = 0.172, and = -0.414, 
which shows that for realistic scattering in a system of 
low symmetry many components of T are nonzero, ex- 
cepting, of course, transport perpendicular to the plane 
of the structure. At the boundary one must study also 
the spin current according to the alternative definition 
14 1 . Previous work [3] has shown this to be of the 
same order of magnitude as the conventional spin cur- 
rent, with occasional sign differences. As our symmetry 
analysis holds for both definitions of the spin current, we 
expect the results to conform to the pattern found above. 

Since several spin components can flow in the same di- 
rection and the same spin component can flow in different 
directions, an interesting experiment could detect spin 
currents. Taking the QW along [113] considered above, 
an electric field applied along x will produce a nonequilib- 
rium spin density, a spin-Hall current and a longitudinal 
spin current composed of spin- a; only. If the QW is joined 
to a material in which no nonequilibrium spin density is 
generated, then a Kerr rotation |3j] or magnetic circu- 
lar dichroism [41| experiment will give a nonzero signal 
exclusively due to the injected spin- a; current. 

In summary, we have shown that in systems with re- 
duced symmetry spin currents are not restricted to the 
spin-Hall effect, and that this fact can help one distin- 
guish experimentally between electrically-induced spin 
densities and spin currents. We have demonstrated in 
addition that spin currents in an electric field are asso- 
ciated with spin precession, whereas spin densities are 
associated with the absence of spin precession. 
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